Entropic Spectral Broadening in Carbon Nanotube Resonators 
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We simulated the behavior of suspended carbon-nanotube (CNT) resonators over a broad range 
of temperatures to address the unexplained spectral broadening and frequency shifts seen in exper- 
iments. We find that thermal fluctuations induce strong coupling between resonance modes. This 
effect leads to spectral fluctuations which readily account for the experimentally observed quality 
factors Q 100 at 300 K. Using a mean field approach to describe entropic fluctuations we analyt- 
ically calculate Q and frequency shifts in tensioned and buckled CNTs and find excellent agreement 
with simulations. 
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Carbon nanotube (CNT) resonators show great 
promise for fundamental science and nanomechanical 
applications owing to their high stiffness, low mass, 
electrical detectability, and defect-free structure. Al- 
though experimental work has shown CNT resonators 
to be highly tunable [1 as well as functional as RF 
transceivers [2 and atomic mass detectors [3 , 4 , these res- 
onators have consistently exhibited much poorer than ex- 
pected properties [U |5]. Specifically, the quality factor 
Q, the key parameter measuring the degree to which an 
oscillating mode is decoupled from its environment, is 
typically less than Q ~ 100 at room temperature. 

These low Qs are not consistent with known dissipa- 
tion mechanisms seen in other nanomechanical resonator 
systems [6 . Analytical phonon-phonon scattering stud- 
ies establish a theoretical upper bound on Q in CNTs 
that is well above experimental values, with Q ^ 50, 000 
at 300 K [7j 8 . Molecular-dynamics simulations of short 
CNTs (1/ ^ 50nm) show interesting behaviors in can- 
tilevered and free CNT segments [9-12 due to anhar- 
monic atomic potentials, and those that make an ex- 
plicit determination of quality factor in thermal equi- 
librium llQi lUj give Q ~ 1,000 at room temperature. 
Typically, however, experiments are performed on CNTs 
with L ^ 1 /im, and the anharmonic elastic effects that 
dominate in short length CNTs do not contribute as sig- 
nificantly at this longer length scale. 

A key difference between CNT resonators and other 
nanomechanical resonators is their nanometer-scale cross 
section [T3^, yielding a very small bending rigidity for flex- 
ural modes and hence large amplitude fluctuations in 
thermal equilibrium. The effects of these thermal fluc- 
tuations on the physics of ID elastic objects has a long 
and storied history. A key length-scale that characterizes 
their properties is the persistence length ip = , where 

is the bending rigidity. Short ID structures {L <^ £p) 
such as microfabricated NEMS are dominated by elas- 
tic effects and thus behave like rigid rods, while long 
ID structures {L ^ ip) such as organic polymers are 
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FIG. 1: a) Schematic of CNT resonator geometry. CNT 
is suspended between two electrodes (in yellow) and has a 
controllable downward force Fz electrostatically applied by a 
voltage to the underlying substrate. A small Fz is applied in 
b) to keep the buckled tube oriented vertically, b) Simulated 
spectral power density of the z-motion of the nanotube as a 
function of strain eo at T = 100 K. The tube is taken from 
tensile strain on the left to compressive strain on the right, as 
illustrated by the insets. Eigenmode shapes are inset, as are 
labels for the eigenfrequencies. 



dominated by configurational entropic effects caused by 
thermal forces and thus take on random 3D configura- 
tions. Micron-scale CNTs behave as semi- flexible poly- 
mers {L < £p)[14 , where the bending energy and con- 
figurational entropy contribute comparably to the total 
free energy. The interplay of these effects in semi-flexible 
polymers is an active area of research. In overdamped, 
aqueous environments, these effects are now well under- 
stood [15 , but high-amplitude fluctuations on the prop- 
erties of nanoscale resonators has not been explored. 

In this Letter, we report simulations of the thermally 
driven dynamics of CNT resonators over temperatures 
from 3 K to 300 K and in the limits of both compressive 
and tensile strain. By simulating purely continuum elas- 
tic behavior, we isolate the specific influence of entropic 
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forces on measured quantities, including both quality fac- 
tor and thermally-induced frequency shifts. Quantum 
mechanical effects are neglected, as the thermal occupa- 
tion of all modes considered are in the classical limit. 
We find that the fluctuations of many resonance modes 
strongly influence the quality factor of CNT resonators, 
an effect we call entropic spectral broadening. We find 
that this entropic spectral broadening is in good agree- 
ment with observed quality factors in experiment. 

CNTs were modeled as ID elastic objects, with bend- 
ing rigidity = tt^^ and extensional rigidity K = irCd, 
where C ^ 345 J/m^ [16] is the 2D elastic modulus of 
graphene and d is the tube diameter. Each CNT was 
discretized into 100 masses joined by axial and torsional 
springs. For a given set of tube dimensions, boundary 
conditions, and externally applied forces, the equilibrium 
geometry was computed via a relaxation method and 
zero-temperature eigenmodes were computed by diago- 
nalizing the force constant matrix [17 . 

In order to study dynamics at non-zero temperature, 
a finite time-difference calculation was employed. Ther- 
mal equilibrium was reached by coupling the nanotube to 
an external heat-bath via a generalized Langevin equa- 
tion, applying a stochastic white noise force with external 
damping. Dynamics at equilibrium were then simulated 
with this Langevin equation using the 4^^ order Runge- 
Kutta method. Adiabatic conditions were also studied 
using Stoermer's rule [18]. The model's validity was con- 
firmed by its ability to properly exhibit equipartition of 
energies for each eigenmode. 

Throughout this Letter, data are shown from simu- 
lations of a typical experimental case of a CNT with 
L = 3 jam and d = 2 nm. We investigate resonance 
properties of this CNT at different zero-temperature 
strains 6o, temperatures T, and externally applied forces 
Fz by quasi-statically sweeping at most one parame- 
ter and measuring the power spectral density of the 
mean ^-displacement (defined in fig la) of the nanotube: 



where tm is the finite 
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time over which the Fourier transform is taken and E 
denotes an ensemble average. The mean z-displacement 
was the chosen parameter to analyze as it is typically 
measured in experiments [19 . From Sz{f)^ we measure 
thermal frequency shifts A/ = f — f(T = OK) and qual- 
ity factor Q = This definition of Q is employed as 
the linewidth Sf is directly measured in experiments. 

We first study the qualitative behavior of this CNT 
as a function of strain at 100 iC. In fig lb, we plot a 
logarithmic color map of Sz (/) with red as the high- 
est amplitude, and blue the lowest. In the tensioned 
limit, to the left of fig lb, the eigenmodes, labeled by 
their modeshape, tune like a tensioned string, giving a 

frequency: f^^^ where ji is the linear mass den- 
sity and N = Ke is the axial tension. Modes without 
mean z-displacements, including the even in-plane (z- 
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FIG. 2: a) Power spectral density of a tensioned nanotube 
with eo = 2 X 10~^. Six temperatures are plotted from 50- 
300 K spaced evenly in T. b) is a linear plot of Sz{f) for 
fiz showing that the frequency shift A/ and linewidth 6f 
both scale linearly with T c) inverse quality factor of fiz vs. 
temperature, d) dimensionless fixed ratio of A/ to Sf that 
corresponds to the number of independent fluctuating modes 
contributing to spectral broadening and frequency shifts. In 
c) and d) data are circles, and theory is a line. The line in c) 
corresponds to Eqn. 2 and in d) rif in Eqn. 1 

direction) eigenmodes and all out-of-plane (x-direction) 
eigenmodes, are not visible. Thus, only modes with 
n = 1, 3, 5... are seen. 

With negative strain, the nanotube undergoes a Eu- 
ler buckling transition, as illustrated in the upper right 
corner of fig lb. With the resulting built-in slack, the 
CNT can bend without stretching leading the spectral 
lines to no longer tune significantly with strain. The 
linewidths, however do tune with strain. Also, there are 
emergent spectral lines corresponding to the motion of 
out-of-plane modes, as labeled to the right of fig. lb. We 
analyze these temperature-dependent linewidths in the 
tensioned case and buckled case separately and discuss 
the different nonlinear coupling mechanisms [^20, ,21^ that 
apply to each. 

We begin with the tensioned case. We select a fixed 
tensile zero-temperature clamping condition eo = 2 x 
10~^, and simulate the nanotube motion at 6 temper- 
atures, incrementing by 50 K from 50 K to 300 K and 
plot in fig. 2a the spectral density over a frequency range 
spanning up to the n = 5 mode. The frequencies of the 
modes shift linearly with T and, as shown in Fig. 2c. The 
peaks also broaden with Q'^ - T. At 300K, Q - 40. 

The quality factor and frequency shifts can be under- 
stood to arise from the change in length of the CNT 
caused by thermal fluctuation in each eigenmode. The 
^th eigenmode is given by Unaiv, t) = an(^)^'na(|;), where 
is the dimensionless modeshape with unit RMS dis- 
placement, an is a time varying amplitude function, and 
both the amplitude a and index a refers to either the 
X OT z direction. The tube elongates by the length 

/^nafj where Pna = SQ^'nai^Y^^ cach independent 
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eigenmode. This leads to the strain being modified as: 

12 \ 

(e) ^ eo + Xln a Pna2r^ where brackets denote a time av- 
erage over the period of oscihation [22 . Assuming that 
{an) obey Boltzmann statistics in thermal equilibrium 
and that fluctuate incoherently, the flrst and second 

moment of the strain shift are Ae = Pna and 
= Y.n,aPla " fl^^ respectively, with thermal 

averages = and = 2 (^^^^ • Here kna is the 
effective spring constant: kna ^ Pna^ + OLna^ where 
ana = /o^^na(^)^^^ parameterizes the bending associ- 
ated with the n*^ eigenmode. In the case of a tensioned 
string, pna ^ n^TT^ and ana ^ n^7r^. 

In the high tension limit, a large number of the res- 
onance modes contribute signiflcantly to length fluctua- 
tions. High frequency modes, however do not contribute 
signiflcantly, as their high bending stiffness limits their 
thermal amplitudes. Approximating the inflnite series as 
an integral gives: 

where the value of Ae is solved self-consistently and n/ 
can be interpreted as the number of independent degrees 
of freedom with signiflcant fluctuation amplitudes. Ae is 
a length-dependent thermal expansion parameter that is 
distinct from the intrinsic thermal expansion simulated 
for CNTs [9\ From Eqn. 1, ~Ke leads to A/ - T . 

Next, we calculate the strain variance: cr^ = ^^sp • 
This leads to the prediction that Sf ~ T. The resulting 
dimensionless quantity 4 ln(2) (^jf^ plotted in flg. 2d is 
a direct measurement of n/. Furthermore, the predicted 
5f constitutes entropic spectral broadening and gives a 
unique prediction for Q: 

in y 

The dominant T-dependence comes from ^ but 
both rif and e have weak T-dependence as a result of 
entropic stretching. The quantity n'jip is the persistence 
length for a tensioned beam and parameterizes how far 
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fluctuations drive the CNT out of equilibrium. The n| 

dependence then can be understood to arise from an ad- 

1 

ditional linewidth broadening ^ Tij due to averaging over 
Uf uncorrelated degrees of freedom. Eqn. 2 accurately 
describes the numerical results, shown in fig. 2c, and 
shows that entropic spectral broadening can account for 
the experimentally observed Q ~ 100 at room tempera- 
ture. 

Next, we study a nanotube under compressive strain, 
picking eo = —4 x 10~^. Applying a small downward 
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FIG. 3: Spectral power density of buckled nanotube with 
eo = — 4 X 10~^ with Fz = 0.6 pN at logarithmically spaced 
temperatures. The built-in intrinsic damping dominates the 
apparent linewidth at the lowest temperature. Inset is a plot 
of of the fsz mode. Simulation data are shown as red 
dots, and theoretical predictions from Eqn. 3 are plotted as 
a blue line. Labels beneath each spectral line correspond to 
their expected origin. 

force Fz = 0.6 pN that mimics the force applied by the 
gate in experiments, simulations were performed at loga- 
rithmically spaced temperatures from 3.1 i^^ to 300 as 
shown in fig. 3. In contrast to the tensioned regime, there 
is complex structure in the spectra, with many spectral 
features growing non-linearly with T. Linear theory pre- 
dicts that fsz and f^z would be the only visible spectral 
lines in fig. 3. The other emergent modes can be identi- 
fied as either oscillations of out-of-plane modes produc- 
ing ^-displacement at twice their natural frequency {2 fix, 
2f2x, and 2fsx) and mixes of in-plane modes {f^z — fsz 
and f^z fsz)- Focusing on the lowest observed in-plane 
mode (fsz), we observe ~ T (fig. 2 inset) and a de- 
crease in frequency with increasing temperature. From 
this we extract Q ~ 5 for this mode at 300 K. 

Further insight is gained by smoothly tuning frequen- 
cies by quasi-statically varying the magnitude of the elec- 
trostatic force Fz = \C'Vg, as shown in fig. 4a, where 
Vg corresponds to the electrostatic voltage applied to the 
gate illustrated in fig lb and C = ^ is the derivative of 
the tube-gate capacitance assuming a 400 nm gap. The 
spectral lines tune with force and the frequency-doubled 
out-of-plane spectral lines cut across the in-plane spectral 
lines. At linear theory predicts there to be no cou- 
pling. However, at the intersections, there are observable 
avoided crossings, even at 10 K as seen in Fig. 4a. This 
behavior is consistent with experiments. To study this, 
we fix Fz at a crossing indicated by the dotted line in fig. 
4a, quasi-statically sweep the temperature, and measure 
Sz{f) (shown in the inset of fig. 4a). Here, we observe 
a sub-linear T-dependence of the splitting at the avoided 
crossing. At 300 K the frequency splitting is nearly 1/5^^ 
of the resonance frequency, indicating a strong coupling 
strength induced entirely by thermal fiuctuations. 
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FIG. 4: a) color-map of spectral noise density at 10 K of the 
mean z-displacement of a buckled nanotube resonator as a 
function of gate voltage assuming tube is 400nm above the 
gate. Modes are labeled as defined in fig. 3. T-dependence 
from K io 300 K of the apparent avoided-crossing at 
150 MHz is shown in the inset at a fixed Vg indicated by 
the red dotted line b) 300 K data of same force sweep in 
a), c) plots the theoretical spectral noise density, based on 
the generalized version of Eqn. 3 described in the text, for 
the same conditions as in b). The dashed lines in b) and 
c) correspond to the theoretical central frequencies from the 
generalized version of Eqn. 3. 

We build an analytical model for this behavior by first 
modeling the CNT as an inextensible object. In this 
framework, when an out-of-plane mode has finite am- 
plitude, the CNT must move up in the z-direction to 
preserve length. To quantify the amount that the tube 
moves in z, we enforce the constraint that the differen- 
tial length dL = Piz \e\dzi — (^ixXn 2ero, 
where the equilibrium deflection of the tube is Zi(^iz{j). 
£^iz is a linear combination of in-plane modes, requiring 

2nd Q^^g^ ^Yi^^ ^^^^ ^ Y PmxXl, whcrC 

_ i^izl^nz) — 2 — which the brackets are an inner 



to 

Vn 



product defined by {a{x)\b{x)) = a{x)b{x)dx [22]. 

Focusing on the interaction of one in-plane and one 
out-of-plane mode, we change the out-of-plane variable 
to Zop = £^op apply time-dependent perturbation 
theory [23] to generate a coupled set of linearized equa- 
tions: 
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dependence, as shown by the theoretical prediction that 
is overlaid on top of the data in the inset of fig. 4a. 

By extending Eqn. 3 to the full interaction matrix, 
populating each mode according to Boltzmann statistics, 
and weighting the resulting frequency probability distri- 
bution by the squared- amplitude distribution we are able 
to generate a theoretical Sz{f) map, which is shown in 
fig. 4c. It compares well with the simulated results at 
300 K shown in fig. 4b. In addition, the linewidths in 
the theoretical Sz{f) give an accurate prediction of the 
simulated Qs, as shown in the inset of fig. 3. The specific 
degree of coupling is geometry dependent, and thus not 
fully analytically generalizable, but Eqn. 3 predicts for 
the lowest in-plane mode fsz that: 
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Eqn. 4 explains the strain dependence of the spectral 
fluctuations seen in the right half of fig. la: Q improves 
with increased buckling as the geometric coupling be- 
tween modes decreases. Furthermore, due to its coupling 
with higher frequency modes, the lowest in-plane mode is 
predicted to decrease frequency with increasing tempera- 
ture at low Vg^ which is frequently observed in experiment 
[1119]. 

These results show that entropic spectral broadening 
dominates the behavior of CNT resonators over a broad 
range of temperatures, and appears to be the main cause 
of temperature-dependent quality factors measured in 
both tensioned and untensioned resonators. To date 
there is limited experimental data characterizing Q over a 
broad temperature range, but thus far, data remain at or 
below our theoretical upper bound [5] fT9] [24 ] . This work 
implies a fundamental limit on Q in high aspect ratio res- 
onators at finite temperatures. Tailoring the geometry of 
these systems is necessary to mitigate thermally-induced 
mode-coupling and thus improve Q. 
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where a{T) = — ^ ^ and (^^p) is time- averaged over 
the period of oscillation. 

Using (xlp) = gives a{T) ^ VT, which when 

substituted into Eqn. 3 predicts average frequency shifts. 
At an avoided crossing, kip ~ Akop gives the prediction 
that A/ ^ ^a{T)fip. This matches the simulated T- 
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